Anisotropic e~^e pressure due to the QED effect in strong magnetic fields and the 
apphcation to the entropy production in neutrino-driven wind 

Kazunori Kohri 

Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto, 606-8502, Japan 

Shoichi Yamada 

Institute of Laser Engineering (ILE), Osaka University, Osaka 565-0871, Japan 



o 
o 

(N 
in 



O 



I 

Oh 



X 



Shigehiro Nagataki 

Department of Physics, School of Science, the University of Tokyo, 7-3-1 Hongo, Bunkyoku, Tokyo 113-0033, Japan 

We study the equation of state of electron in strong magnetic fields which are greater than the critical 
value Be — 4.4 x 10^^ Gauss. We find that such a strong magnetic field induces the anisotropic 
pressure of electron. We apply the result to the neutrino-driven wind in core-collapse supernovae and 
find that it can produce large entropy per baryon, S ~ 400fcB. This mechanism might successfully 
account for the production of the heavy nuclei with mass numbers A = 80 - 250 through the r-process 
nucleosynthesis. 
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Recent years strongly magnetized neutron stars {B ^ 
10^^ Gauss) called "magnetar" have been reported, (see 
recent compilations ||l|,^). Their magnitude is much 
stronger than the critical value Be = m1/e ~ 4.4 x 10^'^ 
Gauss if we assume that their spin-down originates from 
the energy loss of the electromagnetic dipole radiation ^ . 
Such a strong magnetic field is intriguing in the astro- 
physical context because unlike the weak magnetic fields, 
some extraordinary phenomena would be expected from 
large energy splitting of Landau levels, large refractive 
indices of photons, photon splitting effect and so on, 
(see Refs. and references therein). In this letter, we 
study the equation of state of electron in strong magnetic 
fields and show that the electron pressure can become 
anisotropic. So far, another papers have been published 
where they have given an application to the anisotropic 
collapse of the neutron gas HQ. We take a different 
treatment in this study and show that the anisotropic 
pressure can generate entropy and that it can have an 
implication for the r-process nucleosynthesis in the neu- 
trino driven wind in core-collapse supernovae. Such a 
large entropy production is required for the successful r- 
process nucleosynthesis in the wind to account for the 
observational solar abundances 

In quantum electrodynamics (QED), the Lagrangian 
density of electron is given by M 
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C = ip[ijd~ nie] -0, 



(1) 



where is the four-component Dirac spinor of electron, 
nie is electron mass, and the slash means the contraction 
with Dirac's gamma matrices /j? = "Ylnil^PtJ.- Then, the 
energy-momentum tensor density is obtained by, 
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where fj,,v = 0,1,2,3 and the metric is 77^1/ — 
diag(-|-I, — f , — f , — f ). From Euler-Lagrange equation. 



with Eq. (^ we obtain the equation of motion of the 
free electron, i.e., free Dirac equation, [i/) — me]tp — 0- 
In a external magnetic field, we know the following pre- 
scription to obtain the correct Dirac equation, id^ — > 
iDfj^ = idfi — eA^, where e denotes the electromagnetic 
coupling constant, and = {Aq,A) is the vector po- 
tential. Namely we have [i/D — me] "i/) = in the mag- 
netic field. When we choose a vector potential such that 
A° = = = 0, and = Bx, the magnetic field is 
B = (0, 0, B), and we obtain the solution, 
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with (j) = e'^Pyy+P-""^ fnCa and X = • (P - eA)(t)/{E 
TOe), where E is the energy of electron, p — {px,Py,Pz) is 
the momentum of electron, cr^'s are Pauli matrices (j = 1, 
2, 3), n denotes the index of Landau levels (n = 0, 1, • • •), 
and a denotes the spin index (a = 1,-1). The two- 
component spinors C,a are represented by C,i — (J) , and 
C_i = (?). And, /„ = f/v/2"n!V^expKV2]^f„(0, 
where ^ = V eB {x — Py/{eB)), and Hn{£,) is a Hermite 
polynomial function. Then, the electron energy is ex- 
pressed by 



E = ^/ml +pI + eB{2n + 1 - a) 



(4) 



From Eq. (j^, we see that the magnetic field breaks the 
equipartition among p^, (py — eBx)^, and p^, because 
the parallel component to the magnetic field pz can have 
an arbitrary value while the perpendicular components 
have only quantized values as (p^) ~ {{py — eBx)^) = 
eB{n -\- 1/2), where the expectation value of an operator 
is defined by {6) = J 'ip'< O'ipd^ x / J 'ip'f'ipd^x. 

Using the above solution, wc obtain each component of 
the energy-momentum tensor density, T'^ = ip [ij'^Diy] ip. 
Adopting the usual normalization in the quantum field 
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theory, we can integrate T'^ and obtain the energy- 
momentum tensor which is defined by 



(5) 



It is easy to check that all the ofF-diagonal compo- 
nents vanish completely. The diagonal components T^\, 
= dia.g{E,—Px,—Py,~'Pz) give the anisotropic pressure 
which should appear in the hydrodynamic equations. 
They can be read off as 

^ ~ , 1 Q , eB ~ pi 

P. =Py = {n+ ---) — , P^ = -^- (6) 

Hereafter we call them "dynamic pressures" . Compared 
to the case of a weak magnetic-field limit, we easily see 
that only the z-component of the dynamic pressure is 
unchanged. For example, Pz becomes i?/3 in the ultra 
relativistic limit. However, the perpendicular component 
(x-, y-component) of the dynamic pressure takes quan- 
tized values and is different from the one in the case of a 
weak magnetic- field limit. 

To discuss the dynamic pressure averaged in statistical 
mechanics, we consider the grand canonical ensemble and 
introduce the grand potential in a unit volume in a 
strong magnetic filed. 
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where (3 ~ l/lksT), ks is Boltzmann's constant, and fi is 
the chemical potential of electron. The thermodynamic 
pressure and the number density of electron are obtained 
by P = —ft, and N = —d^/dfj,, respectively. Then, we 
find that the equation of state is P = NT even in the 
strong magnetic field. For the thermodynamic pressure 
and the number density of positron, we can obtain them 
only by changing the signature of the chemical poten- 
tial /i — > —fj, in the case of electron. Using Eq. (^), we 
can calculate the dynamic pressure averaged in the grand 
canonical ensemble. 
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where i runs x, y, and z, Si = 2 — 6oi, and Ei = 
^ml +pI + 2eBl. 

Next, we consider the entropy production as a con- 
sequence of the anisotropic dynamic pressure and apply 
it to astrophysical conditions. Here we discuss the r- 
process nucleosynthesis. It has been a mystery for some 
time where and how heavy nuclei whose mass number 
is A = 80 - 250 are produced in the universe. Such 
heavy elements are supposed to be synthesized via so- 
called r-process. The currently favored site is a neutrino 
driven wind in core-collapse supernovae. This is because 



there are a lot of free neutrons near the surface of neu- 
tron star. It has been reported, however, that we need 
a very high entropy per baryon S ^ 400A:b for the suc- 
cessful r-process nucleosynthesis while both numerical 
simulations and analytical treatments fall short of the re- 
quirement by a factor of 2 - 3 . Here we show that we 
can realize such a high entropy per baryon in the strong 
magnetic field by way of anisotropic dynamic pressure, 
and it can be a viable candidate for successful r-process 
nucleosynthesis . 

The first law of thermodynamics is represented by. 
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+ P§r[l/{p/mN)\, where S is the en- 
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tropy per baryon, T is the temperature, mj^ is the nu- 
cleon mass, t is the time, U = em^y is the internal en- 
ergy per baryon (e : specific internal energy), and p 



is the energy density of the matter. = do 
means the Lagrangian time derivative, where Vj is the 
velocity of the fluid (j = 1,2,3). Hereafter the re- 
peated suffix means the summation. P is the thermody- 
namic pressure which is given by the grand potential as 
P = — il. The continuity equation -^p = pdjvj is trans- 
formed into (l/p) = ^/pdjVj. The Euler's equation 
of motion, do{pVi) + dj{pVjVi) = djPij with the conti- 
nuity equation is integrated as (|Dp/2) = VidjPij/ p, 
where P^ is the stress tensor of the dynamic pressure 
(= diag(— Pa;, —Py,—Pz) ) which corresponds to the spa- 
tial parts of the energy-momentum tensor. The equa- 
tion of the energy conservation, do[p{\v\'^/2 + U/rrifq)] + 
di[p{\v\^ /2 + U/mN)vi — PijVj] = is transformed into 
(|i;|V2 + C//mAr) = 1/ pd,{P,jVj). From the above 
equations, we obtain the time-evolution equation of the 
entropy. 



DS ^ P-P, ^ 
T— = ^ — d,v,. 
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Here we emphasize that the pressure in showing up the 
equation of motion is not the thermodynamic pressure 
P but the ensemble averaged stress tensor Pij. This is 
understood as follows. The hydrodynamic equation can 
be obtained from the Boltzmann equation by integrating 
over the momentum. The Boltzmann equation, on the 
other hand, can be obtained from the Green function by 
the so-called gradient expansion as 

df [p' - eA\x)] df e [p' - eA^{x)] dA^{x) df 
dt me dx^ nie dx^ dp^ 

= 0, (10) 

in the current case, where / is the distribution func- 
tion jl^, (see also ]l^ , p^ ). Assuming a local thermal 
equilibrium and integrating Eq. (^) after multiplying 
[p~ eA] , we obtain the magneto- hydrodynamic equations 
with the stress tensor given by Eq. (^). 

To apply the above formulation to the neutrino-driven 
wind in core-collapse supernovae, we adopt a simple an- 
alytic wind solution without a magnetic field which is 



2 



employed in Ref. |T^] where a steady state flow and spher- 
ical symmetry are assumed. The dynamic equations are 
given by 



dv _ IdP GM 

dr p dr r^ ' 

de P dp 

dr p^ dr 



(11) 
(12) 

(13) 



where M is the constant mass-outflow rate in the ejecta, 
r is the radial coordinate from the center of the neu- 
tron star, V is the radial outflow velocity, P is the to- 
tal thermodynamic pressure, G is the Newton's gravita- 
tional constant, and M is the mass of the neutron star. 
The neutrino heating rate of nucleons is represented by 
QeN, where q^N is a neutrino absorption 
rate on free nucleons {v + N ^ e + N'), is a elastic 
scattering rate off background electrons [v + e v + e), 
and (jeN is an electron absorption rate on free nucleons 
{e + N ^ v + N') At the temperature T > 0.5MeV, 
the thermal bath is constituted of 7, e'^e" and nucleons. 
Then, the thermodynamic pressure and the specific inter- 
nal energy are represented by P = IItt'^T'^/ 180 + pT/mN 
and e = Il7r2r''/(60p) + 3T/(2mAr), respectively. On the 
other hand, at the temperature T < O.SMeV we assume 
q—0 because electrons and positrons disappear, and free 
nucleons are bound into a-particles or heavier nuclei. 

To solve the above set of dynamic equations, we should 
give both the boundary and initial conditions. Here we 
assume that the radius of the neutron star Rq is equal to 
the neutrino sphere R^. Then, we give the boundary and 
initial conditions there as M = 1.4Mq, Rq = R^ ~ 10 
km, and p = lO^^g/cm"^. The initial velocity is cho- 
sen so that M does not become more than the crit- 
ical value Merit because the wind should be subsonic, 
e.g., M ^ O{l0-^)Mo erg/s We also assume that 
the neutrino luminosities are identical for all the neu- 
trino species, i.e., Li, = 3.6 x lO^^erg/s and take as a 
mean energy of each neutrino species e^^ = 12 MeV, 
ep^ = 22 MeV, and e^^ = 34 MeV. As for the initial 
temperature, from the analytical treatments we estimate 
Tj ~ 3 MeV for the present parameters |10|. In addition, 
we fix a value of Kg to the final one which is estimated by 
Ye J ~ 0.43 in these parameters 10 1 through the whole 
period in which the radial coordinate evolves from Rq to 
the radius of the outer boundary (~ 10^ km) because the 
dynamics is not sensitive to Ye very much. 

In Fig. H we plot the evolutions of the physical quanti- 
ties as a function of the radial coordinate r. Because we 
are now interested in the entropy, it is sufficient to inves- 
tigate them until r < 100 km which corresponds to the 
temperature T 0.5 MeV. Here we refer to the entropy 
in this model as ^std because this is a standard case in 
a weak magnetic-field limit. 



To investigate the entropy production in strong mag- 
netic fields, we integrate Eq. using the wind solu- 
tion of the above set of dynamic equations. We see from 
Eq. that the wind parallel to B does not generate 
extra entropy. Therefore, we consider a case in which 
the wind flows along the x- or y-axis and use the radial 
coordinate r instead of x or y. Then, we can simplify 
Eq. into one dimensional form and obtain the total 
increment of the entropy in a strong magnetic fleld as 



AS = 



P| - P| dv 
pT/rriN dr 



dt. 



(14) 



where the dynamic electron-positron pressure is deflned 
by the summation, P^ = Pf + PI , which are presented 
in Eq. (|) as Pf ^Pr{p = fie) andPf = P^{^i = -fie)- 
Note that P^ is exactly equal to P*^ which is the ther- 
modynamic pressure appearing in the equation of state, 
that is, P'^ = NeT obtained from the grand potential 
in Eq. (|^), where Ne = Ne- + Ne+- The number den- 
sities of electron and positron should satisfy the follow- 
ing condition of the chemical equilibrium with proton, 
Ne- — Ne+ — Np by which the chemical potential of elec- 
tron /ie is actually determined. In addition, note that 
Eq. ( p^ gives the positive increment of the entropy in 
the wind solution. 
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FIG. 1. Plot of the evolution of the physical quantities as 
a function of the radial coordinate r. They are the entropy 
per baryon number SsTo/ks (long dashed line), the radial 
velocity v in units of 10^ cm/sec (dot-dashed line), the tem- 
perature T in MeV (dashed line), and the energy density of 
the matter p in lO^g/cm'' (solid line). 

Now we compute the cases of strong magnetic fields in 
which V eB 3> T > nie- Then, it is expected that there 
exists large anisotropy of the dynamic pressure. For sim- 
plicity, it would be adequate for us to consider only the 
ground state and the next one in the Landau levels (n = 
0,1) because of the Boltzmann suppression of the Fermi 
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distribution. In Fig. || we plot the evolution of the en- 
tropy in a strong magnetic field. Here we assume that 
the configuration of the magnetic filed is represented by 
B = Bo{r / Ro)~™- , with Bq the magnitude of the mag- 
netic field at the surface of the neutron star (= 5 x lO'^^ 
Gauss) . The index m = 3 corresponds to the dipole mag- 
netic field. From the figure we find that the anisotropic 
dynamic pressure can produce the extra entropy AS/kB 
which is much larger than the standard value. In Fig. ^ 
we plot the extra entropy production as a function of 
the magnetic field at the surface of the neutron star. 
From the figure, we can see that if the magnetic field is 
i? ~ (4 — 6) X 10^^ Gauss, we have a large entropy which 
is required for the successful r-process nucleosynthesis. 
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FIG. 2. Plot of the evolution of the extra increment of the 
entropy in a strong magnetic filed (solid line). For the sake 
of comparison, we also plot the standard case in the weak 
magnetic- field limit S'std/^s (dashed line). 

In this study we have investigated the equation of state 
of electron in strong magnetic fields B ^ Be {— '^A x 10^^ 
Gauss) and found that it induces the anisotropic dynamic 
pressure of electron. We have applied the anisotropic dy- 
namic pressure to the r-process nucleosynthesis in core- 
collapse supernovae. As a result, we obtain a large en- 
tropy S ~ 400A;b for B ~ (4 - 6) x lO^^ Qauss. This 
mechanism of the entropy production might successfully 
give an account of the observational solar abundances of 
the heavy nuclei through the r-process nucleosynthesis 
in the magnetized neutrino wind. Though in this let- 
ter we have assumed for simplicity that the magnetic 
field would not influence the dynamics of the wind, we 
definitely need more consistent magneto-hydrodynamic 
models of the wind. That will be discussed in a separate 
paper p[. 
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FIG. 3. Plot of the extra entropy production AS/ks as a 
function of the magnitude of the magnetic field at the surface 
of the neutron star. The dashed line represents the case of 
m = 2 which corresponds to the configuration which means 
the flux conservation of the magnetic field. The solid line 
represents the case of m = 3 which corresponds to the config- 
uration of the dipole. 
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